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CHANNEL ESTIMATION IN COMMUNICATIONS 

RELATED APPLICATIONS 
[0001] The present Patent Application claims priority to Provisional Patent 
Application Nos. 60/710,372 and 60/797,039, filed 08/22/2005 and 05/01/2006, 
respectively, which are each assigned to the assignee hereof and filed by the 
inventors hereof and which are incorporated by reference herein. 

BACKGROUND 

I. Field 

[0002] The present invention relates generally to telecommunications, and 
more specifically to techniques for channel estimation in a communication 
system. 

II. Background 

[0003] In a typical telecommunications system, the data to be transmitted is 
encoded with a coder, such as, for example, a turbo coder, which generates a 
sequence of symbols, referred to as "code symbols." Several code symbols 
may be blocked together and mapped to a point on a signal constellation, 
thereby generating a sequence of complex "modulation symbols." This 
sequence may be applied to a modulator, which generates a continuous time 
signal, which is transmitted over a wireless channel. 

[0004] At the receiver, the demodulator generates a sequence of modulation 
symbols or soft decisions. If soft decisions are generated, each soft decision 
represents an estimate of a modulation symbol that was transmitted over the 
channel. The estimates may be used to compute the log-likelihood ratio (LLR) 
of the code symbols. In the case of turbo coding, the turbo decoder uses the 
sequence of code symbol LLRs in order to decode the data that was originally 
transmitted. 

[0005] When computing the LLRs of the code symbols, the propagation 
conditions of the channel should be considered. The channel conditions, or the 
channel impulse response, may be estimated at the receiver from a known pilot 
sequence embedded in the data transmission. By way of example, in 
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Orthogonal Frequency Division Multiplexing (OFDM) systems, a Least Squares 
procedure is often used to estimate the channel. Using this procedure, the 
channel may be estimated from a set of pilot tones uniformly spaced across the 
frequency band of one OFDM symbol. 

[0006] However, this channel estimate will be noisy because the pilot tones 
are corrupted with noise. Typically, in order to suppress the noise in the channel 
estimate for OFDM demodulation, the receiver typically averages (i.e. "filters") 
the channel estimate from multiple OFDM symbols. In other words, the channel 
estimate is time-averaged across OFDM symbols. 

[0007] When the channel estimate is time-averaged, the noise in the channel 
estimate is suppressed. In other words, the noise variance of the time-averaged 
channel estimate is smaller than that of the channel estimate without time- 
averaging. The longer the time-averaging window, the more the noise 
suppression. However, the time-averaged channel estimate, even without its 
noise component, may not be the same as the true channel in the presence of 
channel variation. This difference is called the channel estimation bias. Hence, 
the channel estimation error includes of the noise and the bias. 
[0008] E.g., if actual channel variation occurs on a timeframe of about the 
duration of one OFDM symbol and a time-averaging window is about the 
duration of 10 OFDM symbols, then the channel estimation will not be accurate 
for the channel even without the noise. 

[0009] E.g., if the actual channel does not vary for about the duration of 10 
OFDM symbols and the time averaging window is about the duration of 4 OFDM 
symbols, then the channel estimation bias is zero. 

[0010] Hence, the longer time-averaging window leads to smaller noise 
variance, but larger channel estimation bias. The shorter time-averaging window 
leads to smaller channel estimation bias, but larger noise variance. 
[0011] Thus, there are two competing goals: (1) The time-averaging window 
should be long enough to sufficiently suppress the noise, but (2) the time- 
averaging window should be short relative to the actual channel variation 
timeframe to avoid excessive channel estimation bias. The problem is how to 
select an optimum or useful time-averaging window in light of the competing 
goals. 
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[0012] The optimal degree of averaging depends on the SNR and the 
channel Doppler, both of which are unknown a priori. 

[0013] When operating in the low SNR regime, longer time averaging 
reduces noise variance at the expense of increasing estimation bias due to 
channel variation over the averaging window. Since noise is the dominant 
source of degradation at low SNR, it is better to have long time averaging in this 
case. 

[0014] If the noise variance is small compared to the channel estimation 
bias, it may be preferable to shorten the averaging duration thereby reducing 
bias in the channel estimate. 

[0015] The optimal degree of averaging depends on the channel Doppler 
and SNR. Ideally we desire a receiver that automatically varies the degree of 
time averaging, in order to obtain the least overall degradation across all 
channel conditions. 

[0016] There are several adaptive channel estimation algorithms. See, for 
example, "Adaptive channel estimation with velocity estimator for W-CDMA 
receiver," in Proc. IEEE VTC, May 2000, pp. 2024-2028, which is incorporated 
herein by reference. Typically, these algorithms estimate the access terminal 
speed and choose the amount of time-averaging based on a mathematical 
model for the fading statistics. 

[0017] A common mathematical model for the fading statistics is Jakes' 
fading process. However, the many simplifying assumptions behind this model 
may not necessarily apply in the "real world". Estimating access terminal speed 
and choosing the amount of time-averaging based on a highly-simplified 
mathematical model may result in significantly inaccurate estimates of the 
actual rate of channel variation. For example, the access terminal could 
determine that the channel is varying much faster than it is in reality, and hence 
it may use a much smaller degree of time-averaging than necessary. 
Conversely, the access terminal could determine that the channel is varying 
much slower than it is in reality, and hence it may use a much larger degree of 
time-averaging than is desirable. In either situation, the discrepancy can 
produce higher channel estimation error, leading to higher decoding error rates. 
This undesirable property is sometimes referred to as "model mismatch" or 
"model sensitivity". 
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SUMMARY 

[0018] Channel estimation is commonly a difficult problem in communication 
systems. Specifically, selecting a useful or optimal time-averaging window of the 
channel estimates is a major hurdle. 

[0019] A solution for an orthogonal frequency division multiplex (OFDM) 
system is to estimate the channel differently for each tone or tap in the channel. 
This can result in a better channel estimate resulting in higher spectral 
efficiency. 

[0020] Fundamentally, the channel variation rate and SNR is different for 
different frequencies or for different delay times. The following is described in 
terms of the time domain, that is, in terms of delay times (hereinafter, simply, 
delays, or channel tap), but the methods and devices described could similarly 
by performed in the frequency domain for OFDM systems. 
[0021] In a first option, channel estimates are computed with a set of 
candidate time-averaging filters, calculating a metric for each of the candidate 
time-averaging filters, selecting a time-averaging filter based on the metrics 
calculated for the candidate time-averaging filters, and using the channel 
estimate of the selected time-averaging filter. 

[0022] The above selection process can be performed separately for each 
delay, also known as tap, or for each frequency. For example, a time averaging 
filter is selected separately for each delay. 

[0023] In a second option, useful or optimal time averaging filters are 
computed directly for two or more of the several taps. A least square based 
algorithm can be used to calculate the time-averaging filter. The algorithm can 
be performed separately for each delay, also known as tap, or for each 
frequency. For example, a time averaging filter is computed separately for each 
delay. 

[0024] Various aspects and embodiments of the invention are described in 
further detail below. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

[0025] The features and nature of the present invention will become more 
apparent from the detailed description set forth below when taken in conjunction 
with the drawings in which like reference characters identify correspondingly 
throughout and wherein: 

[0026] FIG. 1 shows a block diagram of a cellular wireless communication 
system; 

[0027] FIG. 2 is a graph illustrating a channel estimation process for a zeroth 
OFDM symbol; 

[0028] FIG. 3 is a graph illustrating simulation results for a 10 km/h access 
terminal; 

[0029] FIG. 4 is a graph illustrating simulation results for a 30 km/h access 
terminal; 

[0030] FIG. 5 is a graph illustrating simulation results for a 120 km/h access 
terminal; 

[0031] FIG. 6 is a diagram illustrating an access point and an access 
terminal. 

[0032] FIG. 7 is a flow diagram illustrating a method for calculating a channel 

estimate in a communication system as described herein; 

[0033] FIG. 8 is a flow diagram illustrating a method for calculating a channel 

estimate in a communication system as described herein; 

[0034] FIG. 9 is a diagram illustrating a demodulator for calculating a 

channel estimate in a communication system as described herein; and 

[0035] FIG. 10 is a diagram illustrating a demodulator for calculating a 

channel estimate in a communication system as described herein. 
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DETAILED DESCRIPTION 

[0036] Motivation 

[0037] To improve the quality of channel estimates in OFDM demodulation 
the receiver typically averages the channel estimate from multiple OFDM 
symbols. The optimal degree of averaging depends on the SNR and the 
Doppler, both of which are unknown a priori. Ideally we desire a receiver that 
automatically varies the degree of averaging, in order to obtain the least overall 
degradation across all channel conditions. 
[0038] Channel Estimation without Time-Averaging 
[0039] Referring to FIG. 1, a series 100 of OFDM symbols is shown. Each 
OFDM symbol 101a, 101b, 101c, etc. includes data tones 103 and pilot tones 
105. Let { h ' ls] } denote the channel impulse response (CIR) for the s-th OFDM 
symbol. The pilbt tones from the s-th OFDM symbol are used to obtain a noisy 
estimate of^^: h ils] = h l [s] + Ji,is] ^ where n,[s] js tne estimation error with 
mean zero and variance ° 2 . In other words, h,[s] is a channel estimate of the I- 
tap of the CIR seen by the s-th OFDM symbol, obtained without time-averaging. 

Channel estimates ^ 107a, ^ 107b, ^ 107c, W 107d, 107e 

and 107f are shown, corresponding to OFDM symbol s = -2, -1, 0, 1, 2, 
and 3, respectively. The estimation error is assumed independent across the 
OFDM symbols. 
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[0040] Next, we estimate the channel coefficients at all the tones by a 
discrete Fourier transform (DFT) onf !|U] }. FIG. 2 is a graph illustrating a channel 
estimation process for a zeroth OFDM symbol. We can let 209 denote the 
channel coefficient at the k-th tone. Accordingly, we have 

tf t [j] = ff t [s] + JV t [s] > 

where 

and L is the delay spread 21 1 assumed by the receiver. We assume that L is no 

smaller than the actual delay spread. 

[0041] Channel Estimation with Time-Averaging 

[0042] The receiver averages the channel estimates across multiple OFDM 
symbols to improve the overall quality of the estimates. The averaging can be 

done in either the time-domain on^'^h or in the frequency domain on& M ). 
Here we focus on time-domain averaging. 

[0043] Without loss of generality we focus on the estimation of the channel 
^ of the t-th OFDM symbol. Suppose that we have a bank of P candidate 



output of the p-th candidate filter as: 

KM= I> P (")fy[*+«] 

[0044] Thus, for each candidate filter p, h ^ is obtained by averaging one 

or more values of h t [s]. That is, h ^ is obtained by averaging values of h> ^ 

according to the filter weights w " (n) of the given candidate filter p. The mean- 
square error (MSE) of the channel estimate for the l-th tap and p-th candidate 
filter is: 



0) 
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[0045] Useful or Optimum Time-Averaging Filter 
[0046] It is desirable to choose a filter which produces a channel estimate 
that is "close to" the true underlying channel in a least-squares sense. In 
particular we seek a candidate filter which minimizes (1). However, since the 

true underlying channel is unknown, (1) cannot be computed directly. 
[0047] One solution to this problem is to a) derive an alternate metric that 
can be computed without knowledge of the underlying channel, and b) show 
that a candidate filter which minimizes this metric also minimizes (1). Expand 
the right side of (1) as follows: 

E\m?-^l E faWKMh E \KM 2 (2) 
[0048] Note that the first term is independent of p and the third term can be 
determined from the candidate filter's output. Therefore, if we can compute the 
second term we can find the candidate filter which minimizes (1) based on the 

second and third terms of (2). Still, the second term depends on which is 
unknown. Hence we rewrite the term E ^ h as follows: 

E{Qk\t\ - *,[*])* W'l ) = E{h;[t]h p ,[t] } - w p (0)<t 2 (3) 
[0049] Substituting this in (2), we get 
MSE Pfi = E^tf - 2Re[i?{/? [t)h pJ [t] }]+ E\h pJ [tf+ 2w p (0)<7 2 
(4) 

[0050] If we know the noise variance 0-2 , the MSE can be decomposed into 
an unknown term that is common to all the candidate filters, and terms that can 
be computed from each candidate filter's output without knowledge of the true 
channel. Consider the cost function defined by removing the common term 



E \ h M from (4) and adding the common known term ^P^l 
M t (p) = MSE^+E M| 2 - \h t [tf | 

= E\\j\[t]-h p ,[tf} + 2w p (0)a 2 



(5) 
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[0051] The candidate filter p that minimizes the metric M ' (/?) also 
minimizes MSEpJ . 

[0052] For a practical implementation, the ensemble average in (5) can be 
replaced by the time average, and the noise variance 0-2 can be replaced by an 
estimate &1 derived from the pilot tones. The noise variance may be estimated 
using an estimation procedure similar to co-assigned U.S. Patent Application 
No. 11/047,347, filed on Jan. 28, 2005, entitled "Noise Variance Estimation in 
Wireless Communications for Diversity Combining and Log-Likelihood Scaling," 
which is incorporated herein by reference. The metric at the l-th tap from the p- 
th candidate filter is then approximately 
M,( J p)«M I '(p) + 2w p (0)^ 2 > (6) 
where 



T 

and the overline denotes averaging over an exponential window of length avg . 
When the channel for the t-th OFDM symbol is estimated, we update M,ip) 
using: 

m x (p) Ji- ±- k oo + j- (\h t m m| 2 + 2w p (0)<j 2 ) 

avg J avg ^ 

where & * is the estimate of the noise variance 0-2 available at the time of 
update. In other words, M iW can be calculated using a one-tap auto- 
regressive filter. 

[0053] Note that if we used M ' {p) as a metric, the minimizing filter would be 

degenerate and there would be no averaging, i.e.,^*'^ ~ ^ . The key twist is 
the addition of 2w ' (0)<? \ and this results in a metric that correctly trades off 
reduced bias, between an instantaneous and averaged channel estimate, 
against increased noise variance. 

[0054] Based on the metric M > {p) , we suggest two methods to determine the 
time averaging filter in the following. 
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[0055] Filter Selection Algorithm (FSA) 

[0056] In a filter selection algorithm there is a bank of different candidate 
averaging filters, which may be, for example each tuned for a certain range of 
SNR and Doppler. For the l-th tap of the CIR (0<;<z,-i) ; we choose the 
candidate filter p -« with the smallest metric: 

p scl (0 = arg min M, (q) 

[0057] For example, suppose that there are three candidate filters for FSA. 
(P = 3), with ^ C } =2 ' N(+) - 2 . in the example, the first candidate filter (p=1) is 
defined as follows: w 1 (-2) = Wl (-l) = Wl (0) = Wl (l) = w 1 (2) = l/5 ; tne second 
candidate filter (p=2) is defined as follows: 



w 2 (-l) = w 2 (0) = w 2 (l) = l/3;w 2 (-2) = w 2 (2) = 0. and the thjrd cand i da t e filter 
(p=3) is defined as follows: ^)-^(~^^)^)-^ = ^ . Note that 



the third candidate filter is a zero filter, having coefficients all equal to zero. 
[0058] In the above example, for the zeroth OFDM symbol, the time 
averaged filter outputs are computed as follows: 
fty [0] = 0.2^ [-2] + 0.2ft, [-1] + 0.2/1, [0] + 0.2ft, [1] + 0.2ft, [2]; 



ft 3 ,,[0] = 0. 

[0059] The metric is computed for each filter and for each tap. For each tap, 
the candidate filter with the lowest metric is chosen as the averaging filter for 
channel estimation for that tap. 
[0060] Least-Square Based Algorithm (LSA) 

[0061] In a least squares based algorithm, we directly solve for the filter that 
minimizes M,( - p>> , i.e., 



where ^0) = [h l [t-N M l---JhUl---,h l [t+N w ] r_ 

w is a real vector of the filter weights, and superscript T denotes a matrix 
transpose operation. Define A = N " +1 , and a vector * 4 as a vector of 0's with a 1 
in the A-th position. 

[0062] The optimum filter for the /-th tap (os / si-i) can be approximated by 



M0]= 



0.33ft, [-1] + 0.33/i, [0] + 0.33ft, [1]; and 




WO 2007/024934 
w opt = R, 1 p, 



11 



(7) 



PCT/US2006/032899 



where R ' is an empirical correlation matrix for the l-tap defined by 
r, =Re{h,(t)h;'(t) } and p, js an emp i r j ca i cross-correlation vector for the l-th 
tap defined by P' =Re fe" wfi ' W ) - ^ 2 ^ . When the channel for the t-th OFDM 
symbol is estimated, update R < and p ' by the following auto-regressive steps: 




where &1 is the estimate of the noise variance ^ available at the time of 
update. In other words, R ' is computed using h, (t + * )h" it + s) , for all and 
P' is computed using h.;(t+s)h,(t+s) , for all 

[0063] After each tap of the CIR is estimated by using the filters determined 
by the algorithms (FSA or LSA), the channel estimate at each tone is obtained 
by performing a DFT of these estimated taps. 
[0064] Explanation Regarding PER-TONE ADAPTATION 
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[0065] The following is a practical example demonstrating why using 
different filters for different tones or taps is advantageous. When time-averaging 
is performed in the frequency domain, it is easy to see why a different filter is 
necessary across different tones. Although each multipath component may fade 
at a similar rate, the frequency domain channel coefficients will generally 
experience different rates of change. As a simple example, consider a channel 
with two equal strength paths. Suppose that the channel H(z) is initially given by 
i+z", and a slot later by (i+^+a-fl*" 4 . The response at f = 0 does not change 
at all, whereas the response at the band-edge f = 0.5 changes from 0 to 2S . The 
different rates of change between f=0 compared to f=0.5 would result in larger 
channel estimation error if the same time averaging filter was used at f=0 and 
f=0.5. 

[0066] Explanation Regarding PER-TAP ADAPTATION 
[0067] When time-averaging is performed in the time domain, the optimum 
filter is usually different from tap to tap because each tap in the CIR has 
different energy. In a typical situation, some of the taps within the assumed 
delay spread L may contain virtually no signal energy. In this case it is 
preferable to "zero out" those taps thereby reducing noise and improving overall 
channel estimation quality. A common solution is to force to zero those taps with 
energy below a certain threshold. This approach is referred to as "hard 
thresholding". But if a tap contains some signal energy, zeroing a tap improves 
channel estimation only if the MSE at that tap exceeds the tap's signal energy. 
Since it is difficult to find an optimum level for hard-thresholding that minimizes 
channel estimation error, it is desirable to have an alternative way to zero out 
noisy taps. 

[0068] Explanation Regarding Soft Thresholding 
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[0069] Recall that a filter based on (6) minimizes the MSE at each tap. 
Suppose now that a certain tap contains less signal energy than other taps. The 
chosen filter at this tap will have smaller coefficients but at the same time, in 
order to reduce noise variance, it will have a larger degree of time-averaging 
than at other taps. As a result, this tap will be effectively de-emphasized 
compared to other taps in the CIR. It follows that this algorithm automatically 
reduces the contribution of the noisy taps to the frequency-domain channel 
estimates, leading to lower MSE in the channel coefficients at each tone. In 
summary, with the per-tap adaptation algorithm, a tap is de-emphasized in an 
optimal sense of minimizing its contribution to overall MSE — per-tap adaptation 
achieves optimal "soft-thresholding". 
[0070] Optimum Hard Thresholding 

[0071] Now suppose that we were only able to implement a fixed filter, yet 
we desired performance similar to hard-thresholding. One solution is to consider 
a bank consisting of the fixed filter and another filter with zero weights, and to 
apply per-tap FSA. This algorithm will optimally select between averaging a tap 
with the fixed filter or forcing it to zero. 
[0072] Simulation Results 

[0073] Simulations were produced using a system similar to Enhanced 
BCMCS described in "Updated harmonized Enhanced Broadcast Multicast 
System Proposal for 1xEV-DO", 3GPP2 TSG-C contribution, C30-2005041 8- 
039R1, Apr. 18, 2005, which is incorporated herein by reference. Each OFDM 
symbol includes 320 tones, and every fifth tone is a pilot tone. Hence, there are 
64 uniformly spaced pilot tones. Among the remaining 256 tones, 240 tones are 
data tones and 16 tones are null tones in the guard band. Four OFDM symbols 
make up an "OFDM slot". Since adjacent slots may not be OFDM slots in the 
system under simulation, channel estimation for symbols within a particular 
OFDM slot is performed using only the CIR estimates from the same slot. 
Hence, the two OFDM symbols on the edges of the slot and in the middle of the 
slot use different filters, respectively. 
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[0074] A 1 .8 Mbps data rate (3072 bits/slot, spectral efficiency of 3.2 bps/Hz) 
with 64-QAM and a dual antenna receiver was simulated. The receiver 
truncates the CIR estimate to L = 32 taps. The fading is modeled as Rayleigh 
with a Pedestrian B ITU channel model, as defined in "1xEV-DO Broadcast 
Multicast Service Evaluation Methodology", 3GPP2 TSG-C contribution, C30- 
20041018-015, Oct. 18, 2004, which is incorporated herein by reference. The 
noise estimate a 2 is obtained from the pilot tones and the channel estimates 
{H k [s]} for those pilot tones. 

[0075] Three techniques were simulated and compared. In the first technique 
simulated, an evaluation was made for selecting between different fixed filters 
for different symbols. The fixed filters are as follows: 

[0076] (1) For the zeroth OFDM symbol, N (+) = 3 and N (_) = 0. The set of 
candidate filters is: Filter 1 has weights 0.9, 0.1, 0 and 0. Filter 2 has weights 
2/3, 1/3, 0 and 0. Filter 3 has weights 0.5, 0.25, 0.25 and 0. 
[0077] (2) For the first OFDM symbol, N (+) = 2 and N (0 = 1. The set of 
candidate filters is: Filter 1 has weights 0.25, 0.5, 0.25 and 0. Filter 2 has 
weights 1/3, 1/3, 1/3 and 0. Filter 3 has weights 0.25, 0.25, 0.25 and 0.25. 
[0078] The candidate filters for the second OFDM symbol are the same as 
the candidate filters for the first OFDM symbol, except reversed in time. The 
candidate filters for the third OFDM symbols are the same as for the zeroth 
OFDM symbol, except reversed in time. The above structure of different filter 
sets is used because of the OFDM slot structure of the simulated 
communication system, in which adjacent slots may not be OFDM slots. The 
zeroth and third OFDM symbols are on the edges of slots, whereas the first and 
second OFDM symbols are in the middle of slots. The methods and devices 
described herein can be applied to communication systems where the same 
sets of filters are applied to all OFDM symbols. 

[0079] Roughly speaking, Filter 1 is best at high Dopplers whereas Filter 3 is 
best at low Dopplers. Filter 2 has been used in the BCMCS evaluation of "1xEV- 
DO Broadcast Multicast Service Evaluation Methodology", 3GPP2 TSG-C 
contribution, C30-20041018-015, Oct. 18, 2004, which is incorporated herein by 
reference. 
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[0080] In a second technique simulated, the effect of hard-thresholding was 
studied by augmenting the above three filters with a zero filter (all weights are 
zero). Finally, in a third technique simulated the benefits of soft-thresholding 
was evaluated using LSA. For all three techniques simulated, 16 slot 
exponential averaging windows with updates every slot were used. 
[0081] The performance of these three techniques is shown in FIGs. 3, 4 and 
5 for user speeds of 10, 30, and 120 km/h. As a baseline for comparison, the 
figures show the performance for a receiver that always uses the same fixed 
filter. In terms of best to worst performance for the fixed filters, the ordering at 
10 km/h is Filter 1, 2, and 3. This ordering is reversed at 120 km/h. Thus, no 
filter is optimum across all Dopplers or speeds. On the other hand, the FSA 
results in performance that is at least as good as the best fixed filter at each 
Doppler. Optimal hard-thresholding yields an additional gain of 0.25 to 0.5 dB 
over selection alone. The LSA outperformes hard-thresholding and provides a 
further gain of 0.25 dB or more at high Dopplers. 

[0082] The channel estimation techniques can be implemented in a wireless 
communication network having access points and access terminals. The access 
points may be, for example, base stations in an OFDM cellular system. The 
access terminals may be, for example, wireless handsets such as cell phones, 
or other terminal devices such as a laptop computer with a wireless 
communication card or a built-in transceiver. 
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[0083] FIG. 6 shows a block diagram of an access point 110 and an access 
terminal 120, for implementing channel estimation as discussed above in 
connection with FIGS. 1 to 5. For the reverse link, at terminal 120, a transmit 
(TX) data processor 314 receives traffic data from a data buffer 312, processes 
(e.g., encodes, interleaves, and symbol maps) each data packet based on a 
selected coding and modulation scheme, and provides data symbols. A data 
symbol is a modulation symbol for data, and a pilot symbol is a modulation 
symbol for pilot (which is known a priori). A modulator 316 receives the data 
symbols, pilot symbols, and possibly signaling for the reverse link, performs 
OFDM modulation and/or other processing as specified by the system, and 
provides a stream of output chips. A transmitter unit (TMTR) 318 processes 
(e.g., converts to analog, filters, amplifies, and frequency up converts) the 
output chip stream and generates a modulated signal, which is transmitted from 
an antenna 320. 

[0084] At access point 110, the modulated signals transmitted by terminal 
120 and other terminals in communication with access point 110 are received 
by an antenna 352. A receiver unit (RCVR) 354 processes (e.g., conditions and 
digitizes) the received signal from antenna 352 and provides received samples. 
A demodulator (Demod) 356 processes (e.g., demodulates and detects) the 
received samples and provides detected data symbols, which are noisy 
estimate of the data symbols transmitted by the terminals to access point 1 1 0. 
For the reverse link techniques described herein may be carried out in Demod 
356. A receive (RX) data processor 358 processes (e.g., symbol demaps, 
deinterleaves, and decodes) the detected data symbols for each terminal and 
provides decoded data for that terminal. 
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[0085] For the forward link, at access point 110, traffic data is processed by 
a TX data processor 360 to generate data symbols. A modulator 362 receives 
the data symbols, pilot symbols, and signaling for the forward link, performs 
OFDM modulation and/or other pertinent processing, and provides an output 
chip stream, which is further conditioned by a transmitter unit 364 and 
transmitted from antenna 352. The forward link signaling may comprise power 
control commands generated by a controller 370 for all terminals transmitting on 
the reverse link to base station 110. At terminal 120, the modulated signal 
transmitted by base station 110 is received by antenna 320, conditioned and 
digitized by a receiver unit 322, and processed by a demodulator 324 to obtain 
detected data symbols. For the forward link, techniques described herein may 
be carried out in demodulator 324. An RX data processor 326 processes the 
detected data symbols and provides decoded data for the terminal and the 
forward link signaling. Controller 330 receives the power control commands, 
and controls data transmission and transmits power on the reverse link to 
access point 110. Controllers 330 and 370 direct the operation of terminal 120 
and access point 110, respectively. Memory units 332 and 372 store program 
codes and data used by controllers 330 and 370, respectively, to implement the 
channel estimation as discussed herein. 

[0086] FIG. 7 is a flow diagram illustrating a method for calculating a channel 
estimate in a communication system as described herein. The method starts at 
step 710. In step 710, a first averaging filter is used to estimate a first portion of 
a communication channel. The first averaging filter may be, for example, Filter 
1 , Filter 2 or Filter 3, described above with respect to FIGS. 3 to 5. The first 
portion may be a first tap or a first tone or frequency as described herein. In 
step 720, a second averaging filter is used to estimate a second portion of the 
communication channel. The second averaging filter may be, for example, Filter 
1 , Filter 2, or Filter 3, described above with respect to FIGS. 3 to 5. The second 
portion may be a second tap or a second tone or frequency as described herein. 
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[0087] FIG. 8 is a flow diagram illustrating a method for calculating a channel 
estimate in a communication system as described herein. FIG. 8 is similar to 
FIG. 7, except that three steps 810, 820 and 830 are added before the steps of 
FIG. 7, which correspond to steps 840 and 850 in FIG. 8. The method starts at 
step 810. In step 810, a first metric for a first candidate filter for the first portion 
of the communication channel is calculated. In step 820, a second metric for a 
second candidate filter for the first portion of the communication channel is 
calculated. In step 830, the first averaging filter is selected from the first 
candidate filter and the second candidate filter based at least in part on the first 
metric and the second metric. 

[0088] Step 810 may be include calculating a first function of a first 
difference between a channel estimate of the first portion of the communication 
channel obtained without time averaging and a channel estimate of the first 
portion from the first candidate filter; calculating a first correction factor as a 
second function of a first noise variance estimate and at least one of a plurality 
of weights of the first candidate filter; adding the first correction factor to the first 
function to obtain a third function; and filtering the third function to obtain the 
first metric. The first function may comprise a magnitude squared function. The 
third function may comprise using an auto-regressive filter to obtain the first 
metric. Step 810 may include calculating the first metric based on a zero weight 
filter having only zeros for a plurality of filter weights of the zero weight filter. 
[0089] The communication system may be, to list two examples, an 
orthogonal frequency division multiple access system or a code division multiple 
access system. 

[0090] The method may also include calculating the first averaging filter 
based at least in part on a first past history of the channel estimates of the first 
portion of the communication channel without time averaging and based at least 
in part on a first noise variance; and calculating the second averaging filter 
based at least in part on a second past history of the channel estimates of the 
second portion of the communication channel without time averaging and based 
at least in part on a second noise variance. 
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[0091] Further, the method may include calculating the first averaging filter 
by: calculating a first outer-product of a first vector formed from the first past 
history of the channel estimates of the first portion of the communication 
channel, obtained without time-averaging, corresponding to a plurality of OFDM 
symbols; filtering the first outer-product to obtain an empirical correlation matrix; 
multiplying the first vector by a complex conjugate of the channel estimate of the 
first portion of the communication channel obtained without time-averaging, to 
obtain a second vector; subtracting a noise variance estimate from one of the 
elements of the second vector, to obtain a third vector; filtering the third vector 
to obtain an empirical cross-correlation vector; and multiplying an inverse of the 
empirical correlation matrix by the empirical cross-correlation vector to obtain 
the first averaging filter. 

[0092] The filtering the first outer-product to obtain an empirical correlation 
matrix comprises using an auto-regressive filter to filter the first outer-product to 
obtain an empirical correlation matrix; and wherein, the filtering the third vector 
comprises using a second auto-regressive filter to filter the third vector to obtain 
an empirical corss-correlation vector. 
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[0093] For the first portion, the first metric may be M,(1) , and the second 
metric may be , as described herein. For the second portion, the first 
metric may be M ^ , and the second metric may be M ^ , as described 
herein. In step 830, the first averaging filter is selected based at least in part on 
the first metric and the second metric of the first portion, and the second 
averaging filter is selected based at least in part on the first metric and the 
second metric of the second portion. For example, the first and the second 
metric of the first portion may be compared and the candidate filter producing 
the lowest metric is selected as a first averaging filter in step 830. Many more 
candidate filters and many more portions are likely in a practical system. Two 
candidate filters, two portions and two metrics per portion are used here for 
illustration. The techniques and systems described will repeated as many times 
as required or convenient for the communication system used, such as and 
OFDM system. In some cases, calculating the first metric for the channel 
estimate of the first portion may include calculating the first metric based on a 
zero weight filter having only zeros for a plurality of filter weights of the zero 
weight filter. 

[0094] In some embodiments, the method includes calculating the first 
averaging filter based on a least squares algorithm (LSA) and calculating the 
second averaging filter based on the least squares algorithm (LSA). 
[0095] FIG. 9 is a diagram illustrating a demodulator for calculating a 
channel estimate in a communication system as described herein. For example, 
FIG. 9 may be included in Demod 356 or demodulator 324, shown with respect 
to FIG. 6. Returning to FIG. 9, module 910 is a module for using a first 
averaging filter to estimate a first portion of a communication channel. The first 
averaging filter may be, for example, Filter 1, Filter 2, or Filter 3, described 
above with respect to FIGS. 3 to 5. The first portion may be a first tap or a first 
tone or frequency as described herein. In module 920, a second averaging filter 
is used to estimate a second portion of the communication channel. The second 
averaging filter may be, for example, Filter 1, Filter 2, or Filter 3, described 
above with respect to FIGS. 3 to 5. The second portion may be a second tap or 
a second tone or frequency as described herein. Modules 910 and 920 may be 
combined into a single module in Demod 356 or demodulator 324. 
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[0096] FIG. 10 is a diagram illustrating a demodulator for calculating a 
channel estimate in a communication system as described herein. FIG. 10 is 
similar to FIG. 9, except that three modules 1010, 1020 and 1030 are added to 
the modules of FIG. 9. For example, FIG. 9 may be included in Demod 356 or 
demodulator 324, shown with respect to FIG. 6. Returning to FIG. 10, in module 
1010, a first metric for a first candidate filter for the first portion of the 
communication channel is calculated. In module 1020, a second metric for a 
second candidate filter for the first portion of the communication channel is 
calculated. For the first portion, the first metric may be M,(l) , and the second 
metric may be M,(2) , as described herein. For the second portion, the first 
metric may be M 2 (l), and the second metric may be M 2 (2), as described 
herein. In module 1030, the first averaging filter is selected from the first 
candidate filter and the second candidate filter based at least in part on the first 
metric and the second metric. 

[0097] For example, the first and the second metric of the first portion may 
be compared and the candidate filter producing the lowest metric is selected in 
module 830. Many more candidate filters and many more portions are likely in a 
practical system. Two candidate filters, two portions and two metrics per portion 
are used here for illustration. The techniques and systems described will be 
repeated as many times as required or convenient for the communication 
system used, such as an OFDM system. In some cases, calculating the first 
metric for the channel estimate of the first portion of the channel may include 
calculating the first metric based on a zero weight filter having only zeros for a 
plurality of filter weights of the zero weight filter. 

[0098] In some embodiments, the method includes calculating the first 
averaging filter based on a least squares algorithm (LSA) and calculating the 
second averaging filter based on the least squares algorithm (LSA). 
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[0099] The techniques and modules described herein may be implemented 
by various means. For example, these techniques may be implemented in 
hardware, software, or a combination thereof. For a hardware implementation, 
the processing units within an access point or an access terminal may be 
implemented within one or more application specific integrated circuits (ASICs), 
digital signal processors (DSPs), digital signal processing devices (DSPDs), 
programmable logic devices (PLDs), field programmable gate arrays (FPGAs), 
processors, controllers, micro-controllers, microprocessors, other electronic 
units designed to perform the functions described herein, or a combination 
thereof. 

[00100] For a software implementation, the techniques described herein may 
be implemented with modules (e.g., procedures, functions, and so on) that 
perform the functions described herein. The software codes may be stored in 
memory units and executed by processors or demodulators. The memory unit 
may be implemented within the processor or external to the processor, in which 
case it can be communicatively coupled to the processor via various means. 
[00101] The previous description of the disclosed embodiments is provided to 
enable any person skilled in the art to make or use the features, functions, 
operations, and embodiments disclosed herein. Various modifications to these 
embodiments may be readily apparent to those skilled in the art, and the 
generic principles defined herein may be applied to other embodiments without 
departing from their spirit or scope. Thus, the present disclosure is not intended 
to be limited to the embodiments shown herein but is to be accorded the widest 
scope consistent with the principles and novel features disclosed herein. 
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CLAIMS 

1. A method for calculating a channel estimate in a communication 
system, the method comprising: 

using a first averaging filter to estimate a first portion of a 
communication channel; and 

using a second averaging filter to estimate a second portion of the 
communication channel, wherein the second averaging filter is different 
from the first averaging filter. 

2. The method of claim 1 , wherein: 

the first portion of the communication channel corresponds to a first 
tap of an impulse response of the communication channel; and 

the second portion of the communication channel corresponds to a 
second tap of the impulse response of the communication channel. 

3. The method of claim 1 , wherein: 

the first portion of the communication channel corresponds to a first 
tone of a frequency response of the communication channel; and 

the second portion of the communication channel corresponds to a 
second tone of the frequency response of the communication channel. 

4. The method of claim 1 , further comprising: 

calculating a first metric for a first candidate filter for the first portion 
of the communication channel; 

calculating a second metric for a second candidate filter for the first 
portion of the communication channel; and 

selecting the first averaging filter from the first candidate filter and 
the second candidate filter based at least in part on the first metric and the 
second metric. 

5. The method of claim 4, further comprising: 
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calculating a third metric for the first candidate filter for the second 
portion of the communication channel; 

calculating a fourth metric for the second candidate filter for the 
second portion of the communication channel; and 

selecting the second averaging filter from the first candidate filter 
and the second candidate filter based at least in part on the third metric 
and the fourth metric. 

6. The method of claim 4, wherein calculating the first metric 
comprises: 

calculating a first function of a first difference between a channel 
estimate of the first portion of the communication channel obtained without 
time averaging and a channel estimate of the first portion from the first 
candidate filter; 

calculating a first correction factor as a second function of a first 
noise variance estimate and at least one of a plurality of weights of the 
first candidate filter; 

adding the first correction factor to the first function to obtain a third 
function; and 

filtering the third function to obtain the first metric. 

7. The method of claim 6, wherein the first function comprises a 
magnitude squared function. 

8. The method of claim 6, wherein the filtering the third function 
comprises using an auto-regressive filter to obtain the first metric. 

9. The method of claim 4, wherein calculating the first metric for the 
first candidate filter for the first portion of the communication channel 
comprises: 

calculating the first metric based on a zero weight filter having only 
zeros for a plurality of filter weights of the zero weight filter. 
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1 0. The method of claim 1 , wherein the communication system is an 
orthogonal frequency division multiplexing communication system. 

1 1 . The method of claim 1 , wherein the communication system is a 
code division multiple access system. 

1 2. The method of claim 1 , further comprising: 

calculating the first averaging filter based at least in part on a first 
past history of the channel estimates of the first portion of the 
communication channel without time averaging and based at least in part 
on a first noise variance; and 

calculating the second averaging filter based at least in part on a 
second past history of the channel estimates of the second portion of the 
communication channel without time averaging and based at least in part 
on a second noise variance. 

13. The method of claim 12, wherein the calculating the first averaging 
filter comprises: 

calculating a first outer-product of a first vector formed from the first 
past history of the channel estimates of the first portion of the 
communication channel, obtained without time-averaging, corresponding 
to a plurality of OFDM symbols; 

filtering the first outer-product to obtain an empirical correlation 

matrix; 

multiplying the first vector by a complex conjugate of the channel 
estimate of the first portion of the communication channel obtained without 
time-averaging, to obtain a second vector; 

subtracting a noise variance estimate from one of the elements of 
the second vector, to obtain a third vector; 

filtering the third vector to obtain an empirical cross-correlation 
vector; and 

multiplying an inverse of the empirical correlation matrix by the 
empirical cross-correlation vector to obtain the first averaging filter. 
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14. The method of claim 13, wherein: 

the filtering the first outer-product to obtain an empirical correlation 
matrix comprises using an auto-regressive filter to filter the first outer- 
product to obtain an empirical correlation matrix; and wherein, 

the filtering the third vector comprises using a second auto- 
regressive filter to filter the third vector to obtain an empirical cross- 
correlation vector. 

15. A device for calculating a channel estimate in a communication 
system, the device comprising: 

means for using a first averaging filter to estimate a first portion of a 
communication channel; and 

means for using a second averaging filter to estimate a second 
portion of the communication channel, wherein the second averaging filter 
is different from the first averaging filter. 

1 6. The device of claim 1 5, wherein: 

the first portion of the communication channel corresponds to a first 
tap of an impulse response of the communication channel; and 

the second portion of the communication channel corresponds to a 
second tap of the impulse response of the communication channel. 

1 7. The device of claim 1 5, wherein: 

the first portion of the communication channel corresponds to a first 
tone of a frequency response of the communication channel; and 

the second portion of the communication channel corresponds to a 
second tone of the frequency response of the communication channel. 

1 8. The device of claim 1 5, further comprising: 

means for calculating a first metric for a first candidate filter for the 
first portion of the communication channel; 

means for calculating a second metric for a second candidate filter 
for the first portion of the communication channel; and 
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means for selecting the first averaging filter from the first candidate 
filter and the second candidate filter based at least in part on the first 
metric and the second metric. 

1 9. The device of claim 1 8, further comprising: 

means for calculating a third metric for the first candidate filter for 
the second portion of the communication channel; 

means for calculating a fourth metric for the second candidate filter 
for the second portion of the communication channel; and 

means for selecting the second averaging filter from the first 
candidate filter and the second candidate filter based at least in part on 
the third metric and the fourth metric. 

20. The device of claim 1 8, wherein the means for calculating the first 
metric comprises: 

means for calculating a first function of a first difference between a 
channel estimate of the first portion of the communication channel 
obtained without time averaging and a channel estimate of the first portion 
from the first candidate filter; 

means for calculating a first correction factor as a second function 
of a first noise variance estimate and at least one of a plurality of weights 
of the first candidate filter; 

means for adding the first correction factor to the first function to 
obtain a third function; and 

means for filtering the third function to obtain the first metric. 

21 . The device of claim 20, wherein the first function comprises a 
magnitude squared function. 

22. The device of claim 20, wherein the means for filtering the third 
function comprises means for using an auto-regressive filter to obtain the 
first metric. 
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23. The device of claim 1 8, wherein means for calculating the first 
metric for the first candidate filter for the first portion of the communication 
channel comprises: 

means for calculating the first metric based on a zero weight filter 
having only zeros for a plurality of filter weights of the zero weight filter. 

24. The device of claim 1 5, wherein the communication system is an 
orthogonal frequency division multiplexing communication system. 

25. The device of claim 15, wherein the communication system is a 
code division multiple access system. 

26. The device of claim 1 5, further comprising: 

means for calculating the first averaging filter based at least in part 
on a first past history of the channel estimates of the first portion of the 
communication channel without time averaging and based at least in part 
on a first noise variance; and 

means for calculating the second averaging filter based at least in 
part on a second past history of the channel estimates of the second 
portion of the communication channel without time averaging and based at 
least in part on a second noise variance. 

27. The device of claim 26, wherein the means for calculating the first 
averaging filter comprises: 

means for calculating a first outer-product of a first vector formed 
from the first past history of the channel estimates of the first portion of the 
communication channel, obtained without time-averaging, corresponding 
to a plurality of OFDM symbols; 

means for filtering the first outer-product to obtain an empirical 
correlation matrix; 

means for multiplying the first vector by a complex conjugate of the 
channel estimate of the first portion of the communication channel 
obtained without time-averaging, to obtain a second vector; 
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means for subtracting a noise variance estimate from one of the 
elements of the second vector, to obtain a third vector; 

means for filtering the third vector to obtain an empirical cross- 
correlation vector; and 

means for multiplying an inverse of the empirical correlation matrix 
by the empirical cross-correlation vector to obtain the first averaging filter. 

28. The device of claim 27, wherein: 

the means for filtering the first outer-product to obtain an empirical 
correlation matrix comprises means for using an auto-regressive filter to 
filter the first outer-product to obtain an empirical correlation matrix; and 
wherein, 

the means for filtering the third vector comprises means for using a 
second auto-regressive filter to filter the third vector to obtain an empirical 
cross-correlation vector. 

29. A communication device comprising: 

a receiver configured to a receive wireless communication signal; 

and 

a demodulator coupled to the receiver and configured to: 

use a first averaging filter to estimate a first portion of a 

communication channel; and 

use a second averaging filter to estimate a second portion of 

the communication channel, wherein the second averaging filter is 

different from the first averaging filter. 

30. The communication device of claim 29, wherein: 

the first portion of the communication channel corresponds to a first 
tap of the communication channel; and 

the second portion of the communication channel corresponds to a 
second tap of the communication channel. 



31 . The communication device of claim 29, wherein: 
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the first portion of the communication channel corresponds to a first 
tone of the communication channel; and 

the second portion of the communication channel corresponds to a 
second tone of the communication channel. 

32. The communication device of claim 29, wherein the demodulator is 
further configured to: 

calculate a first metric for a first candidate filter for the first portion 
of the communication channel; 

calculate a second metric for a second candidate filter for the first 
portion of the communication channel; and 

select the first averaging filter from the first candidate filter and the 
second candidate filter based at least in part on the first metric and the 
second metric. 

33. The communication device of claim 32, wherein the demodulator is 
further configured to: 

calculate a third metric for the first candidate filter for the second 
portion of the communication channel; 

calculate a fourth metric for the second candidate filter for the 
second portion of the communication channel; and 

select the second averaging filter from the first candidate filter and 
the second candidate filter based at least in part on the third metric and 
the fourth metric. 

34. The communication device of claim 32, wherein calculating the first 
metric comprises: 

calculating a first function of a first difference between a channel 
estimate of the first portion of the communication channel obtained without 
time averaging and a channel estimate of the first portion from the first 
candidate filter; 

calculating a first correction factor as a second function of a first 
noise variance estimate and at least one of a plurality of weights of the 
first candidate filter; 
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adding the first correction factor to the first function to obtain a third 
function; and 

filtering the third function to obtain the first metric. 

35. The communication device of claim 34, wherein the first function 
comprises a magnitude squared function. 

36. The communication device of claim 34, wherein the filtering the 
third function comprises using an auto-regressive filter to obtain the first 
metric. 

37. The communication device of claim 32, wherein calculating the first 
metric for the first candidate filter for the first portion of the communication 
channel comprises: 

calculating the first metric based on a zero weight filter having only 
zeros for a plurality of filter weights of the zero weight filter. 

38. The communication device of claim 29, wherein the demodulator is 
further configured to: 

calculate the first averaging filter based at least in part on a first 
past history of the channel estimates of the first portion of the 
communication channel without time averaging and based at least in part 
on a first noise variance; and 

calculate the second averaging filter based at least in part on a 
second past history of the channel estimates of the second portion of the 
communication channel without time averaging and based at least in part 
on a second noise variance. 

39. The communication device of claim 38, wherein the calculating the 
first averaging filter comprises: 

calculating a first outer-product of a first vector formed from the first 
past history of the channel estimates of the first portion of the 
communication channel, obtained without time-averaging, corresponding 
to a plurality of OFDM symbols; 
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filtering the first outer-product to obtain an empirical correlation 

matrix; 

multiplying the first vector by a complex conjugate of the channel 
estimate of the first portion of the communication channel obtained without 
time-averaging, to obtain a second vector; 

subtracting a noise variance estimate from one of the elements of 
the second vector, to obtain a third vector; 

filtering the third vector to obtain an empirical cross-correlation 
vector; and 

multiplying an inverse of the empirical correlation matrix by the 
empirical cross-correlation vector to obtain the first averaging filter. 

40. A machine readable medium having instructions for calculating a 
channel estimate in a communication system stored thereon, the 
instructions comprising: 

a first instruction for using a first averaging filter to estimate a first 
portion of a communication channel; and 

a second instruction for using a second averaging filter to estimate 
a second portion of the communication channel, wherein the second 
averaging filter is different from the first averaging filter. 

41 . The machine readable medium of claim 40, wherein: 

the first portion of the communication channel corresponds to a first 
tap of an impulse response of the communication channel; and 

the second portion of the communication channel corresponds to a 
second tap of the impulse response of the communication channel. 

42. The machine readable medium of claim 40, wherein: 

the first portion of the communication channel corresponds to a first 
tone of a frequency response of the communication channel; and 

the second portion of the communication channel corresponds to a 
second tone of the frequency response of the communication channel. 
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43. The machine readable medium of claim 40, the instructions further 
comprising: 

a third instruction for calculating a first metric for a first candidate 
filter for the first portion of the communication channel; 

a fourth instruction for calculating a second metric for a second 
candidate filter for the first portion of the communication channel; and 

a fifth instruction for selecting the first averaging filter from the first 
candidate filter and the second candidate filter based at least in part on 
the first metric and the second metric. 

44. The machine readable medium of claim 43, further comprising: 

a sixth instruction for calculating a third metric for the first candidate 
filter for the second portion of the communication channel; 

a seventh instruction for calculating a fourth metric for the second 
candidate filter for the second portion of the communication channel; and 

an eighth instruction for selecting the second averaging filter from 
the first candidate filter and the second candidate filter based at least in 
part on the third metric and the fourth metric. 

45. The machine readable medium of claim 43, wherein the third 
instruction for calculating the first metric comprises: 

a ninth instruction for calculating a first function of a first difference 
between a channel estimate of the first portion of the communication 
channel obtained without time averaging and a channel estimate of the 
first portion from the first candidate filter; 

a tenth instruction for calculating a first correction factor as a 
second function of a first noise variance estimate and at least one of a 
plurality of weights of the first candidate filter; 

an eleventh instruction for adding the first correction factor to the 
first function to obtain a third function; and 

a twelfth instruction for filtering the third function to obtain the first 

metric. 
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46. The machine readable medium of claim 43, wherein the first 
function comprises a magnitude squared function. 

47. The machine readable medium of claim 43, wherein the filtering the 
third function comprises using an auto-regressive filter to obtain the first 
metric. 

48. The machine readable medium of claim 43, wherein the third 
instruction for calculating the first metric for the first candidate filter for the 
first portion of the communication channel comprises: 

a thirteenth instruction for calculating the first metric based on a 
zero weight filter having only zeros for a plurality of filter weights of the 
zero weight filter. 

49. The machine readable medium of claim 40, the instructions further 
comprising: 

a fourteenth instruction for calculating the first averaging filter 
based at least in part on a first past history of the channel estimates of the 
first portion of the communication channel without time averaging and 
based at least in part on a first noise variance; and 

a fifteenth instruction for calculating the second averaging filter 
based at least in part on a second past history of the channel estimates of 
the second portion of the communication channel without time averaging 
and based at least in part on a second noise variance. 

50. The machine readable medium of claim 49, wherein the third 
instruction for calculating the first averaging filter comprises: 

a sixteenth instruction for calculating a first outer-product of a first 
vector formed from the first past history of the channel estimates of the 
first portion of the communication channel, obtained without time- 
averaging, corresponding to a plurality of OFDM symbols; 

a seventeenth instruction for filtering the first outer-product to obtain 
an empirical correlation matrix; 
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an eighteenth instruction for multiplying the first vector by a 
complex conjugate of the channel estimate of the first portion of the 
communication channel obtained without time-averaging, to obtain a 
second vector; 

a nineteenth instruction for subtracting a noise variance estimate 
from one of the elements of the second vector, to obtain a third vector; 

a twentieth instruction for filtering the third vector to obtain an 
empirical cross-correlation vector; and 

a twenty-first instruction for multiplying an inverse of the empirical 
correlation matrix by the empirical cross-correlation vector to obtain the 
first averaging filter. 

51 . The machine readable medium of claim 50, wherein: 

the filtering the first outer-product to obtain an empirical correlation 
matrix comprises using an auto-regressive filter to filter the first outer- 
product to obtain an empirical correlation matrix; and wherein, 

the filtering the third vector comprises using a second auto- 
regressive filter to filter the third vector to obtain an empirical cross- 
correlation vector. 
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